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HEAT-TRANSFER ANALYSIS FOR LIQUID-METAL FLOW IN RECTANGULAR 
CHANNELS WITH HEAT SOURCES IN THE FLUID 
by Robert M. Inman 
Lewis Research Center 

SUMMARY , 

For an understanding of turbulent liquid-metal heat transfer in rectangular channels 
with heat generation in the fluid stream, an analysis was performed for forced- convection 
heat transfer to slug flow in rectangular channels with aspect ratios from 1 to °° and with 
uniform heat sources in the fluid. 

The analysis was carried out for both uniform wall heat flux aild uniform wall temper- 
ature. The results obtained apply in the thermal entrance region as well as in the fully 
developed region. For the wall-heat-flux situation, the four channel walls were consid- 
ered to be uniformly heated - the heat flux on the short sides was considered an arbitrary 
fraction or multiple of the heat flux on the long sides. For the uniform temperature case, 
the wall temperature was assumed to be constant and equal on the four walls of the 
channel. 

The analysis was based on the additional assumptions that (1) heat transport by eddy 
conduction is negligible compared with that by molecular conduction, (2) internal heat 
generation is spatially uniform, and (3) fluid properties are invariant with temperature. 
The temperature distributions were determined by utilizing the method of superposition. 

The effects on temperature distributions of (1) the ratios that determine the role of 
internal heat generation relative to that of wall heat transfer or temperature driving 
force, (2) specified aspect ratio, and (3) specified wall heat fluxes or wall temperatures 
around the channel periphery were investigated. Numerical results for wall tempera- 
tures, heat-transfer rates, and bulk mean fluid temperatures are presented graphically 
as functions of significant parameters. The solutions identify the locations of maximum 
temperature or wall heat-transfer rates. 

The results are useful for estimating local heat-transfer characteristics in turbulent 
heat- generating liquid-metal flow in rectangular channels when the Peclet number is 
approximately 100. 



INTRODUCTION 


Technological developments in space power generation have stimulated interest in the 
problem of forced-convection heat transfer to liquid-metal flow in passages with internal 
heat generation in the fluid stream. This system applies, for example, in the design of 
liquid-metal-fuel reactors, electromagnetic pumps and flowmeters, and liquid-metal 
magnetohydrodynamic (MHD) generators. The fluid in these devices will be heated by 
radioactive fission products or by an electric current flowing through the fluid. A sche- 
matic diagram of a simple MHD generator, which consists of a rectangular channel with 
sidewall electrodes, is shown in figure 1. The liquid metal flows in the channel with ve- 
locity U in the x-direction. A uniform magnetic field H is applied in the y-direction. 
Electric current flows from the right electrode through an external load to the left elec- 
trode, and then through the liquid metal back to the right electrode to complete the circuit. 
In many cases, the flow in such devices will be accompanied by heat, both that dissipated 
internally through viscous or Joule heating and that occurring through wall or electrode 
heating. 

For the proper operation of such devices, a satisfactory temperature distribution 
along the passage walls must be maintained. The designer, therefore, must be able to 
compute the temperature distribution along the walls and to know how much heat must be 
removed to cool the walls and thus prevent temperatures from exceeding design limits. 

The problem involves studying liquid-metal channel flow with combined internal heat 
sources and wall heat transfer. Turbulent flows, which are very often encountered in 
practice, are studied in this analysis. 

This investigation is concerned with hydrodynamically developed turbulent flow of a 
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liquid metal in rectangular channels with aspect ratios from 1 to °° and with uniform in- 
ternal heat generation. An analysis including both uniform wall heat flux and uniform wall 
temperature is carried out for both the thermal entrance region and the fully developed 
region. For the case of uniform wall heat flux, the four channel walls are considered to 
be uniformly heated, but the heat flux on the short sides of the channel is an arbitrary 
fraction or multiple of the heat flux on the long sides. The walls of an MHD generator 
(fig. 1), for example, must insulate electrically, and if the wall temperature is less than 
that of the liquid metal, the wall must conduct heat from the fluid into an adjacent coolant. 
In addition, heat produced by electric currents in the electrodes will be transferred to the 
liquid metal. These effects may lead to unequal heat addition on adjacent sides. Unequal 
heat addition may also be the result of unwanted heat leakage or addition through insula- 
tion. For the uniform temperature case, the wall temperature is assumed constant and 
equal on all four walls of the channel. 

This report is concerned with the determination of the axial and peripheral tempera- 
ture distribution and heat-transfer characteristics in rectangular channels with aspect 
ratios from 1 to 00 . Attention is focused on the rectangular channel because of its in- 
creasing use in the applications mentioned. Within the knowledge of the author, the ex- 
perimental and analytical studies of turbulent liquid- metal channel flows with internal 
heat generation have been confined to elementary geometries such as the circular tube 
(refs. 1 to 4) and the parallel- plate channel (refs. 4 to 6). In contrast to this moderate 
amount of information, turbulent-flow heat transfer to a heat-generating liquid metal in 
a rectangular channel has apparently received little analytical and no experimental work. 

A few studies related to the problem considered in the present investigation are 
noted. In the absence of internal heat generation in the fluid, solutions have been devel- 
oped which approximate situations that might occur with liquid metals in turbulent flow 
through rectangular channels. Fully developed slug- flow Nusselt numbers and wall- 
temperature distributions are presented in references 7 to 9. In these references, the 
problem was solved by assuming, in addition to a uniform velocity throughout the channel, 
that turbulent eddying does not contribute to conduction of heat within the fluid. The re- 
sults pertain to systems characterized by low Reynolds and Prandtl numbers and to the 
portion of the rectangular channel beyond the thermal entrance region. In the discussion 
section of reference 10, Hoagland presents work done on the thermal entrance region 
for laminar slug flow in rectangular channels and gives some numerical results. Slug- 
flow wall-temperature distributions in the entrance region of channels having a variety of 
noncircular shapes, including a square channel, are presented in reference 11. The data 
were obtained by an analog technique. In reference 12, forced- convection heat transfer 
to laminar slug flow in a rectangular channel is examined for the boundary condition of a 
channel wall temperature both peripherally and axially uniform. This analysis was 
carried out without consideration of internal heat sources. 
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Magnetohydrodynamic heat transfer in a rectangular channel with specified heat- flux 
boundary conditions is considered in reference 13. The flow is assumed to be hydromag- 
netically and thermally fully developed, and heat conduction is by molecular conduction 
only. The effects of viscous dissipation and Joule heating are included in a general 
manner. 

The present, very limited knowledge of turbulent liquid-metal flow and eddy diffu- 
sivity variation in noncircular passages makes theoretical progress very unlikely without 
appeal to simplified models. It does not seem likely, moveover, that any single model, 
for which a mathematical analysis is feasible, will prove adequate for all Reynolds and 
Prandtl numbers. Therefore, in order to gain some understanding of the complex prob- 
iem of turbulent liquid-metal flow in rectangular channels with wall heat transfer and in- 
ternal heat sources, consideration is given to a simplified but representative model. 

This specific model not only retains many of the physical characteristics of turbulent 
liquid-metal channel flow, but also leads to a tractable mathematical problem. This 
model should, therefore, provide information on the temperature distribution and heat- 
transfer characteristics for such flows in rectangular passages. 

The idealized system assumed to approximate the forced- convection system under 
consideration is based on the following postulates: 

(1) The established turbulent velocity profile is represented by a uniform distribution. 

(2) The thermal eddy diffusivity is small compared to the thermal molecular diffu- 
sivity and is neglected. 

(3) Longitudinal heat conduction is small compared to longitudinal convection and 
transverse conduction and is neglected. 

(4) The internal heat generation is spatially uniform. 

It is pointed out (e. g. , refs. 7 and 14) that the blunt-nosed turbulent velocity distribution 
for a liquid-metal system can be represented satisfactorily by a uniform distribution. 
(References 15 and 16 remark that this assumption is even better satisfied for turbulent 
flow in MHD devices. ) The second postulate implies (ref. 17) that the thermal solution 
pertains to systems characterized by Peclet numbers equal to, or less than, approxi- 
mately 100. The third postulate (ref. 18) introduces a negligible error for Peclet moduli 
that are equal to, or greater than, approximately 100. Turbulent-flow heat transfer to 
liquid metals can be estimated, at least in the absence of internal heat generation, by the 
use of the slug- flow solutions for molecular conduction (refs. 7 to 9). 

In this analysis, numerical results are provided for the case of internal heat sources 
which are uniform across the channel cross section and along the channel length. The 
results can be extended, however, to include sources which vary in the transverse and 
longitudinal directions (refs. 19 and 20). 

This investigation is divided into two major sections dealing, respectively, with in- 
ternal heat generation in the presence of prescribed wall-heat fluxes and with internal 
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heat generation in the presence of an axially and peripherally uniform wall temperature. 
Beginning the analysis with the case of internal heat generation and prescribed wall heat 
fluxes will facilitate the development for the case of internal heat generation with uniform 
equal channel surface temperatures. 


INTERNAL HEAT GENERATION WITH PRESCRIBED WALL HEAT FLUXES 

The rectangular channel and its coordinate system for the present case are shown in 
figure 2. The turbulent velocity is postulated to be fully established at x = 0. The fluid 
temperature at the channel entrance is uniform across the section at a value t g . Within 
the channel a heating process takes place that includes a uniform heat generation within 
the liquid metal and a uniform heat transfer at the channel walls. The fluid is assumed 
to have constant physical properties, and only steady- state heat transfer is investigated. 

The established turbulent velocity profile is represented by 

u = U = constant (1) 

(All symbols are defined in appendix A. ) The differential equation describing convective 
heat transfer for the idealized system takes the form 

U — = ... tL /d. .1 + 1±\ + _Q_ (2) 

9x P c p\dy 2 3z / P c p 



^Uniformly heated 

sidewalls, q s * aq B CD-8776 


Figure 2. - Coordinate system for rectangular channel with different uniform heating on 
each pai r of opposite wall s . 


5 



and the boundary conditions are 


at = _^b 

3y k 


= - — at y = 0 (specified wall heat flux) 


(3a) 


at 

— = — at y = a (specified wall heat flux) 

By k 

9t = .^s 
Bz K 


at z = 0 (specified wall heat flux) 


(3b) 


(3c) 


at _ n s 

3z K 


- — at z = b (specified wall heat flux) 


t(0, y, z) = t (entrance condition) 


(3d) 

(3e) 


It can be shown (appendix B) that the following equation satisfies the differential equa- 
tion (2) and the previous conditions, and is thus the solution to the problem where internal 
heat generation and wall heat transfer occur simultaneously: 


t-te 

q R a 


•[■(■*;) 


+ R 


x' + y' 2 - y' + — z' 2 - ofz' + — (1 + era) 
a 6 


OO 

X 

k=2, 4, . . 


a i_2 2 , 

— — cos(k7ry')e" K n x - aa 

, 2 2 
k ir 


OO 

y ^ 

/ j ,2 2 

t—J l IT 


COS 


(a 5V i2 ’ 2x 'A’ 2 


(4) 


1—2, 4, . 


where 


R = ^ 

«B 


(5) 


There may arise practical applications where it is desirable to express the heat flux 
on the broad walls in terms of sidewall heat transfer as q^ = o;q s . The temperature dis- 
tribution can also be expressed in terms of sidewall heat transfer by substituting 1/a 
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for a, aq s for q B , and replacing the parameter R by the parameter S, which is de- 
fined as 


S = ^ ( 6 ) 

% 

The parameters R and S are the ratios of internal heat evolution to the heat transferred 
at the channel walls, and they give a measure of the importance, in connection with tem- 
perature development, of internal heat generation relative to wall heat transfer. 


Heat-Transfer Results 


From the temperature distribution given by equation (4), various quantities of engi- 
neering interest can be determined. In the developments to follow, the analytical results 
are expressed in terms of the broad-wall heat transfer q B> 

Of practical application are the wall-temperature variations corresponding to pre- 
scribed wall heat transfer and internal heat generation. The local temperatures 
t w (x’, z') = t w B along the broad walls are found from equation (4) by evaluating it at 
y' = 0 or at y' = 1: 


V, B ~ *e 

q B a 


2 1 + 


6 


+ R 


x’ + — z' 2 - az* + A (1 + ora) 
c r 6 


- V e- k2A ’ 
kV 

k=2, 4, . . . 

The local temperatures t w (x*,y’) - t w g along the sidewalls are found in a similar man- 
ner from equation (4) by evaluating it at z' = 0 or at z’ = or: 


OO 
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t - t 
w, s e 


q B a 


= 2^1 + — ^ + rJx* + y' 2 - y* + — (1 + aa) 
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A 2 zZ 
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( 8 ) 


1-2, 4, . .. 


Another more convenient form of these equations is obtained by introducing the bulk 
mean temperature t^. For a uniform heat source and a uniform wall heat transfer, the 
bulk temperature is given by 


*b ~ *e 

q B a 



K 


( 9 ) 


Then the local wall- to bulk-temperature difference along the broad walls is given by 


V B ~ tb = £ z > 2 - az » + l(l + 


(1 + acr) - 


q B a 


OO 
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4 e -k 2 jr 2 x* 


kV 


k=2, 4, 


- aa 


f -J- C os(ln*)e - l2 ” 2 *'/° 2 

zZ A 2 v °> 


( 10 ) 


1 = 2 , 4 ,... 


while the temperature difference along the sidewalls is 
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T 2 -. 2 

6 7T 


(ii) 


Equations (10) and (11) are independent of the internal heat generation rate and, therefore, 
are valid for all heat generation rates for the proposed model. The differences between 
the fully developed wall and bulk temperatures are 


^w, B ~ Wj 9 ! 

1 = “ z' 2 - az' +1(1 + 

q B a o 


(1 + ao) 


^w, s ” 9 , 

q = y' Z - y' + 1 (1 + ao) 

q B a 6 


(12a) 


(12b) 


The ratios of local to fully developed temperature differences along the broad walls 
or sidewalls at any location in the channel are found from equations (10) to (12) as 

OO OO 

-±-cos(u*y i2 ' ,2x ' / ° 2 

u w, B ~ V J _ k=2, 4, . . . Z=2,4, ... 
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(13b) 


It can be observed from the foregoing equations that a prescribed heat flow from the chan- 
hel surfaces to the fluid produces, at a given axial location, a local wall temperature or 
a wall- to bulk-temperature difference which varies around the periphery of the channel 
and which also varies with channel aspect ratio. 


Results and Discussion 

To illustrate the effects of internal heat generation and prescribed wall heat transfer 
on wall-temperature distributions, local bulk mean temperature, and thermal entry 
lengths for liquid-metal flow in rectangular passages, a number of solutions have been 
obtained for various combinations of aspect ratio cr, wall heat flux ratio a or a, and 
heat flux ratio R or S. For brevity, only some of these results are included in this 
report. In particular, values of the parameters a and a chosen for the computations 
correspond to the following cases: 

(1) For a = - 1, the broad walls are heated and the sidewalls are cooled. 

(2) For a = 0, the sidewalls are insulated. 

(3) For a = 1, uniform heating (or cooling) takes place all around the channel 
periphery. 

(4) For a = 0, the broad walls are insulated. 

Designers of channels for space-power-generation systems are interested in the tem- 
peratures reached by the walls, and in particular the peripheral location where the wall 
temperature will assume its highest value for a known wall heat input and a known internal 
heat generation rate. From an examination of the analytical results, it is expected that 
the peak temperatures will occur either in the corner of a channel or at the centerline of 
the broad wall or short wall, depending on the wall heat flux ratio a. Therefore, knowl- 
edge of temperature conditions in such regions is essential. The wall temperatures 
given here are relative to the bulk temperature, since the heat flux ratio R or S is 
then eliminated as a parameter. 

The longitudinal variations of the dimensionless wall temperatures were evaluated 
from the analytical solution, and are presented in figure 3 for channel aspect ratios of 
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Dimensionless wa 





1, 4, and 10. For the special case of insulated sidewalls (a = 0), the solution is independ- 
ent of the channel aspect ratio; therefore, the result applies for all aspect ratios. The 
result corresponds to that result for heat- generating slug flow between parallel plates with 
uniform wall heat flux. At some places along a wall, the wall- to bulk-temperature dif- 
ference may be negative. This is understandable if it is recalled that t is a local value 
along the wall, while t^ is an average value over the entire cross section. When the 
broad walls are heated and the sidewalls are cooled ( a = -1), the peak temperatures occur 
along the centerline of the broad walls at all axial positions for the aspect ratios shown 
(fig. 3(a)). When only the broad walls are heated (a = 0) (fig. 3(b)), these walls assume 
a uniform temperature, for all aspect ratios, which is higher than the temperatures along 
the insulated sidewalls. For a rectangular channel with uniform heat flux around the en- 
tire periphery (fig. 3(c)), the peak temperatures for all cr occur at the channel corners. 
Finally, when heat is transferred from only the short walls (fig. 3(d)), these walls attain 
a uniform temperature that is higher than the temperatures along the adiabatic broad 
walls. 

Another quantity with practical application is the bulk-mean-temperature variation 
along the length of the channel. The bulk temperature is given here relative to the tem- 
perature of the fluid at the entrance to the channel. The dimensionless bulk temperatures 
are presented in figure 4 for channel aspect ratios of 1, 20, and °° (parallel- plate channel) 
with the heat flux ratio appearing as a family parameter. 

Positive and negative values of the heat flux ratios R and S are considered in fig- 
ure 4. It is assumed that Q is positive (a heat source). A positive value of R, there- 
fore, implies that q^ is positive - that is, that heat is being transferred from the broad 
walls to the fluid. A negative value of R, on the other hand, implies that q B is nega- 
tive, or that heat is being transferred from the fluid to the broad walls. For positive R, 
therefore, internal heat generation and broad-wall heat transfer reinforce one another to 
produce a bulk temperature larger than that obtained in the absence of internal heat gen- 
eration. Conversely, for negative R, heat transfer to the broad wall opposes internal 
heat generation in the bulk-temperature development. Similar arguments apply in connec- 
tion with sidewall heat transfer q and heat flux ratio S. 

For very small values of |r| or |s}, wall heat transfer dominates the bulk- 
temperature development, while for large values of [r| or |s[, the effects of internal 
heat generation dominate. This explains the variety of trends that are evident in each of 
the figures. The bulk-temperature development is, in general, affected slightly by chan- 
nel aspect ratios a to about 20 and insignificantly thereafter. For a channel with insu- 
lated sidewalls, however, the bulk-temperature development is independent of the channel 
aspect ratio, and corresponds to that for heat- generating slug flow between parallel plates 
with uniform wall heat flux. 

The foregoing presentation of results has been concerned with wall temperatures in 
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(a) Uniform heating of broad walls and uni- (b) Uniform heating of broad walls (q s - 0). 

form cooling of sidewalls (q s /q B - -1). Channel aspect ratio, 1 < b/a <». 



0 .1 .2 .3 .4 0 .1 .2 .3 .4 


Dimensionless longitudinal distance, (x/a)/RePr 

(c) Uniform heating on four walls (q s ■ q B ). (d) Uniform heating of sidewalls (q B » 0). 

Figure 4. - Longitudinal variation of bulk mean temperature for various values of heat flux ratio Oa/q B 
or Qa/q s . 

the thermal entrance region. As a matter of general interest, the wall-temperature vari- 
ation around the periphery of the channel in the fully developed region is considered. The 
wall temperatures are again given relative to the bulk temperature, since in the fully de- 
veloped region the temperature difference (t - t, ) is independent of x. In addition, as 

W D (J 

noted earlier, the temperature difference t - is independent of the heat flux ratio R 
at all axial positions. 

Fully developed wall temperatures are presented in figure 5 for the values of a and 
~a considered earlier and for various aspect ratios. It is worthwhile to recall that t 
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is a local value along the wall, while is an average value over the entire cross sec- 
tion. Therefore, at some places along the wall, the temperature difference (t - t.) is 

d 

negative, which means that is larger than t . The hot spots are strikingly displayed 
in these figures, appearing in the corners, or at the broad-wall or short-wall midpoints. 

It is also apparent that the aspect ratio has a significant effect on the temperature distri- 
bution. 

Of considerable practical importance to the designer is a knowledge of the conditions 

under which entrance effects must be accounted for in heat-transfer calculations. In 

theory, the approach of a local- to bulk-temperature difference to the fully developed 

value is asymptotic. Consequently, it is difficult to identify a specific length of channel 

as a thermal entrance length. It is practice to define a thermal entrance length in terms 

of the downstream distance x/aRePr at which the temperature difference comes within 

5 percent of the fully developed value. The variation with dimensionless axial distance of 

local to fully developed wall temperatures at the channel corner and at the broad-wall and 

short-wall centerlines was evaluated from the analytical expressions. The results thus 

obtained are plotted in figure 6 for values of cn = - 1 and a = 0 and for parametric aspect 

ratio values. The corner temperature ratio is not shown for the square channel when 

a = -1, because for this situation the corner temperature difference is zero at all axial 

positions. Lines delineating the condition (t - t. )/ (t - t^) =0. 95 have been drawn in 

d 

figure 6 to graphically illustrate the approach to fully developed conditions. 

From an inspection of the graphs, it is seen that, for the wall heating conditions rep- 
resented, the wall-temperature profiles approach fully developed conditions at widely dif- 
ferent distances from the entrance to the heating section, depending on the particular wall 
location chosen for consideration. It is also evident that the sidewalls of the channels 
strongly influence the temperature developments, with the thermal entrance lengths de- 
creased with decreased aspect ratio. 

It is the practice in reporting heat transfer connected with flow through noncircular 
passages to present average heat-transfer coefficients or Nusselt numbers based on a 
heat flow averaged around the channel periphery and on an average wall temperature. 

This practice might be useful when the wall temperature remains constant everywhere, 
or at least is constant around the periphery of a channel at a given axial position. In the 
present situation where the wall boundary condition is one of peripherally and axially uni- 
form heat input, it is apparent that the knowledge of the resulting local wall-temperature 
distribution is more important to the designer than Nusselt numbers or even average 
thermal- entry lengths. Therefore, no attempt has been made to determine the latter 
quantities. It is felt that the results are presented in a form more convenient for engi- 
neering calculations. 
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Wall-temperature ratio, 



(a) Along channel corner with uniform heating (b) Along wall centerlines with uniform heating 
of broad walls and uniform cooling of side- of broad walls and uniform cooling of sidewalls 


walls (q s /qg = -1). tc ls^B = 



0 .1 .2 .3 A 0 .1 .2 .3 .4 

Dimensionless longitudinal distance, (x/a)/RePr 


(c) Along channel corner with heat transferred (d) Along wall centerlines with heat trans- 
from short walls (q B = 0). ferred from short walls (q B = 0). 

Figure 6. - Wall-temperature ratios in thermal entrance region. 
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INTERNAL HEAT GENERATION WITH UNIFORM WALL TEMPERATURE 


Consideration is now given to the problem where a liquid metal with a uniform tem- 
perature t enters a rectangular channel whose four walls are maintained at uniform 
equal temperatures t , which are different from the entrance value (fig. 7). In addition, 



Figure 7. - Coordinate system for rectangular channel with uniform and equal wall temperatures. 


a uniform heat generation begins within the fluid at x = 0. Although many wall- 
temperature combinations are possible, the boundary conditions of axially and peripher- 
cally uniform wall temperature are of some technical interest and will serve to demon- 
strate the effects of internal heat generation on heat-transfer characteristics in the pre- 
sence of uniform wall temperatures. 

The velocity profile is again represented by equation (1), and the energy transfer pro- 
cesses in the liquid metal are assumed to be still governed by equation (2). The boundary 
conditions are now 

t = t at y = 0 and y = a (specified wall temperatures) (14a) 

t = t w at z = 0 and z = b (specified wall temperature) (14b) 

t(0, y, z) = t (entrance condition) (14c) 

A solution to equation (2) that satisfies the boundary conditions of equations (14) is given 

by 
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where 


N = 



(16) 


and m = 1, 3, 5, . . . and n = 1, 3, 5, . . . . The parameter N is essentially the 
ratio of temperature difference t g - t to internal heat generation. The derivation and 
discussion of equation (15) is presented in appendix C. 


Heat-Transfer Results 

From the temperature distribution given by equation (15), quantities of engineering 
interest can be determined. The local bulk temperature T b for the situation of combined 
internal heat generation and uniform equal wall temperatures, for example, is obtained 
from (appendix C) 
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When the channel wall temperatures are specified, the quantity of greatest practical 
interest is the longitudinal variation of the wall heat fluxes that is required to maintain 
the wall temperatures constant. If, for example, the wall temperatures are less than 
that of the liquid metal entering the channel, and in addition internal heat sources are pre- 
sent, the walls must conduct heat from the liquid metal into a surrounding coolant. An 
expression for the local peripheral heat flux along the sidewall z’ = 0 is obtained by 
applying Fourier’s law 




z=0 


w' 


= -«(, -t )[il 


3z 
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to the temperature solution (eq. (15)) to obtain 
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The local heat flux q (cr) along the opposite sidewall z' = cr is obtained as the negative 
of that along the wall z’ = 0; that is, q z (a) = -q z (0). The local peripheral heat flux along 
the broad channel wall y' = 0 is found in a similar manner with the result 
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The local dimensionless heat flux q y (l) at the opposite broad wall y’ = 1 is given by 

q y (i) = -q y (o). 

It should be noted from equations (18) and (19) that the boundary condition of a con- 
stant peripheral and axial channel wall temperature requires (at a given axial position) 
the local heat flow from the channel surfaces to the fluid, or from the fluid to the channel 
surfaces, to vary around the periphery. The local heat flow also varies with channel as- 
pect ratio. 





Dimensionless longitudinal distance, (x/al/RePr 

(c) Channel aspect ratio, 10. (d) Channel aspect ratio, ® (parallel-plate 

channel). 

Figure 8. - Bulk or mixed-mean temperature variation for internal heat generation with equal wall 
temperatures. 
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Results and Discussion 

To illustrate the results, local bulk mean temperatures and wall heat fluxes were 
computed for several combinations of aspect ratio a and the heat generation param- 
eter N. Bulk temperature results for internal heat generation with peripherally and 
axially uniform wall temperatures are presented in figure 8 for channel aspect ratios of 
1, 4, 10, and °° and for parametric values of the heat generation ratio. It is again as- 
sumed that Q is positive, but clearly T g - T w may be either positive or negative, 
depending on whether the walls are cooler or warmer than the entering liquid metal. 

The heat generation ratio or parameter N, therefore, is also either positive or nega- 
tive. Consequently, both positive and negative values are considered in the figure. The 
curves of (t. - t )/ (t - t ) change continuously from an initial value of one to an asymp- 

D W G W oo oo 

tote (for X — °°) of (1/N) EE 64/m^r? n^T^ nn 

m=l n=l 

as t^ varies from t g to 

OO OO 

‘w + £‘e - l w>/ N ] E E 64 / m2 " 2 ' 4l tn- In 1116 

m=l n=l 

absence of internal heat sources (N — °°), the bulk 
temperature asymptotically approaches the temper- 
ature of the walls. 

To illustrate more clearly the influence of chan- 
nel aspect ratio, the dimensionless bulk temperature 
is plotted against the dimensionless longitudinal dis- 
tance in figure 9 for values of 1/N equal to - 10, 0, 
and 10, with the channel aspect ratio as the family 
parameter. The bulk temperature development for 
a given value of 1/N is affected by channel aspect 
ratio cr to about 10 and insignificantly affected 
thereafter. 

A plot of the longitudinal variation of the dimen- 
sionless sidewall heat flux at the wall locations y* = 0. 25 and 0. 50 was prepared and is 
presented in figure 10. Channel aspect ratios of 1, 4, 10, and 20 are considered, and 
the heat generation ratio is the parameter. The heat fluxes commence with large negative 
values near the channel entrance and change continuously with increasing downstream dis- 
tance. For a channel of given aspect ratio, the magnitude and sign of the parameter N 
have a profound influence on the wall heat fluxes. In the absence of internal heat genera- 
tion the wall heat fluxes decrease to zero at large downstream distances. For negative 
values of N, q^ (0) can change sign at some position along the wall length, which indicates 



Dimensionless longitudinal distance, 


Figure 9. - Bulk temperature dependence on 
aspect ratio for various values of heat 
generation ratio. 
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Dimensionless longitudinal distance, (x/a)/RePr 
(c) Channel aspect ratio, 10. (d) Channel aspect ratio, 20. 

Figure 10. - Local sidewall heat-flux results for internal heat generation with equal wall tempera- 





a change in direction of the heat transfer at that location along the wall. For a given 
value of 1/N, the influence of channel aspect ratio on wall heat fluxes is slight for the 
range of aspect ratios chosen for illustration. When a — °° (infinite parallel plates), 
however, the sidewall heat fluxes diminish to zero at all positions along the sidewalls 
(see eq. (18)). 

Longitudinal variations of the dimensionless broad-wall heat fluxes at the wall loca- 
tions z”(= z ’/or) of 0. 25 and 0. 50 were evaluated and are presented in figure 11 for chan- 
nels with aspect ratios of 1, 4, 10, and 20. The heat generation ratio N appears again 
as the parameter. The figure shows the influence of aspect ratio on wall heat fluxes. 

As <7 is increased, for a specified value of 1/N, the peripheral wall heat flux distribu- 
tion becomes more uniform. The broad-wall heat fluxes, at a given longitudinal position 
and for a selected value of N, remain almost constant for all aspect ratios a beyond 
approximately 10. The heat generation ratio has a significant influence on the broad-wall 
heat fluxes for a channel of given aspect ratio. For N — °°, the wall heat fluxes decrease 
to zero at large downstream distances. For negative values of N, it is observed that 
q y (0) also may change sign at some position along the broad- wall length. 

It is the practice in considering heat transfer in noncircular passages to define a 
heat-transfer coefficient as the average heat flux at the wall (averaged around the channel 
periphery) divided by the difference between the fluid bulk temperature and the wall tem- 
perature. In the foregoing figures, however, for negative values of the parameter N, 
the wall- to bulk-temperature difference became zero when the bulk temperature of the 
fluid was equal to the wall temperature and negative when the wall temperature was less 
than the bulk temperature. In addition, for negative values of N, there occurred changes 
in direction of the wall heat transfer at certain locations. Infinities and zeros, therefore, 
occurred in the heat-transfer coefficient, and this variation opens to question the utility 
of a heat-transfer coefficient under the condition of internal heat generation in the fluid. 
Provision of adequate, adjacent cooling channels for effective heat removal to avoid over- 
heating the fluid and to maintain the required wall temperatures is of practical concern 
to the designer. Also, the local wall-heat-flux conditions are necessarily more impor- 
tant from the standpoint of effective heat removal than consideration of the average heat- 
transfer conditions. It is conjectured that independent consideration of local wall- to 
bulk-temperature differences and wall heat fluxes will prove more useful for preliminary 
heat-transfer design purposes. For these reasons, information on the average heat- 
transfer coefficients for rectangular channels containing flowing, heat- generating liquid 
metals are not reported. 
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Dimensionless broad wait heat flux, qy(0)a/K{t e - t w ) 
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CONCLUDING REMARKS 

Theoretical results have been presented for heat transfer to liquid metals flowing in 
rectangular channels with prescribed wall heat fluxes or wall temperatures and heat 
sources uniformly distributed in the fluid. The various boundary and internal heating ef- 
fects considered herein have yielded a qualitative understanding of the axial and periph- 
eral temperature or heat-flux distributions as well as general heat-transfer character- 
istics. 

It should be emphasized that the present results are strictly valid only for the slug- 
flow velocity distribution with heat transfer by molecular conduction. The simplification 
provided by these assumptions, however, has made it possible to obtain exact mathemat- 
ical solutions to the governing energy equation. A fact of greater importance is that 
these temperature distributions approximate conditions expected for turbulent liquid- 
metal flow in rectangular channels for relatively low Prandtl and Reynolds numbers. 

Some of the characteristics of turbulent, heat-generating liquid-metal heat transfer 
are the following: 

1. Wall temperature distributions for specified wall heat fluxes are, in general, 
sensitive to the channel aspect ratio. 

2. For the slug-flow velocity distribution, wall- to bulk-temperature differences for 
specified wall heat fluxes are independent of the internal heat generation rate. It is con- 
jectured, therefore, that with nonuniformity of the velocity profile the influence of an in- 
ternal source on wall- to bulk-temperature differences will be slight. 

3. The thermal entrance length is quite sensitive to channel aspect ratio and wall 
heating conditions. For the conditions investigated, the thermal entrance length is least 
for a square channel. 

4. Bulk temperature development in channels with uniform temperature along all 
walls is influenced by channel aspect ratio to an aspect ratio of about 10. Beyond this 
value, the influence of aspect ratio is slight. 

5. The wall- heat- flux variations required to maintain the wall temperatures constant 
are dependent on the channel aspect ratio, internal heat source strength, and magnitude 
of the entrance temperature in relation to the value of the wall temperature. In partic- 
ular, for wall temperatures cooler than the fluid entrance temperature, a change in di- 
rection of the heat transfer may occur at some position along a given wall. An engineer- 
ing design of a wall heating or cooling system would necessarily have to take this into 
consideration. The results, in addition, point out the locations of maximum wall temper- 
atures or maximum wall heat-transfer rates. 

More refined analyses would undoubtedly take into account velocity- profile and eddy- 
diffusivity variations to provide a more realistic description. Turbulent velocity profiles 
for circular-pipe and parallel-plate channel systems have been conveniently and satisfac- 
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torily represented by power-law expressions. Reference 21 presents and applies a vari- 
ational method for determining velocity distributions for flow of a power- law fluid in rec- 
tangular channels. The results might prove useful for heat-transfer studies. Within the 
knowledge of the author, eddy-diffusivity distributions for turbulent flow in rectangular 
channels have received little theoretical consideration. In reference 6, the writer derived 
analytical heat-transfer solutions for turbulent flow of a heat- generating liquid metal in a 
parallel- plate channel for high Reynolds numbers where the effect of radial eddy- 
diffusivity variations must be included. Uniform wall heat flux and uniform wall- 
temperature cases were considered. A linear approximation of the complicated eddy- 
diffusivity function was used as in references 22 and 23, and good agreement of the heat- 
transfer results with existing numerical solutions was demonstrated. A linear diffusivity 
function could be applied to the problem of turbulent heat transfer in rectangular passages. 
However, within the author's knowledge, no work has thus far been done along these lines. 

For fully developed turbulent flow in straight channels of rectangular cross section, 
there exists a transverse mean flow, commonly known as secondary flow, superimposed 
upon the axial mean flow (refs. 24 to 27). This secondary flow interacts with the axial 
mean flow and the turbulence structure in a complex manner. Qualitative observations in- 
dicate that secondary flow tends to equalize the velocities and temperatures within any 
cross section of a rectangular channel. However, because little quantitative knowledge 
exists on the influence of the secondary flow on temperature distributions within a fluid, 
particularly for fully developed turbulent liquid-metal flow, an exact theoretical calcula- 
tion of the temperature field in rectangular channels is not possible at the present time. 
Consequently, the wall -temperature and wall -heat -flux distributions determined in the 
present analyses may be taken as limiting cases. 

The present analytical treatments are adequate for preliminary design purposes and 
should give the designer a reasonable estimate of the locations and magnitudes of maxi- 
mum temperatures or wall heat fluxes that could be expected for a given channel config- 
uration in the presence of internal heat generation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 1, 1966, 

129-01-09-07-22. 
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APPENDIX A 


SYMBOLS 


An coefficients in series for temper- 

ature distribution with internal 
heat generation and zero wall 
temperatures 

a length of short side 

b length of broad side 

C coefficients in double series for 

inn 

temperature distribution with 
internal heat generation and 
zero temperatures 

Cp specific heat of fluid at constant 
pressure 

D „ coefficients in double series for 
rs 

temperature distribution with 
internal heat generation and 
zero wall temperatures 

H uniform magnetic field 

5 <i * - v 

N heat generation ratio, 

Q 

Pr Prandtl number, pCp//c 

Q rate of internal heat generation 

per unit volume 

q heat flux per unit wall area 

Qy(0) peripheral heat flux along broad 
wall y' = 0 

q z (0) peripheral heat flux along side- 

wall z' = 0 

R heat flux ratio, Qa/q B 

Re Reynolds number, pUa/ju 


S heat flux ratio, Qa/q„ 

D 

T dimensionless temperature, 

< l - - l w> 


t temperature 

U velocity, uniform over channel 



cross section 

u 

local fluid velocity 

W n 

function of y', defined by eq. (C14) 

X 

coordinate measured along the 


axial direction 

X T 

dimensionless axial coordinate, 


x/aRePr 

y 

coordinate measured along short 


side 

y’ 

dimensionless coordinate, y/a 

z 

coordinate measured along long 


side 

z* 

dimensionless coordinate, z/a 

z" 

dimensionless coordinate, 


z/b = z'/cr 

a 

wall heat flux ratio, q s Alg 

a 

wall heat flux ratio, q B /q s = l/ a 

^n 

eigenvalue, nw^l + (l/o^) 

r 2 

mn 

function of m, n, and a, 
rr 2 |m 2 + (n 2 /a 2 )] 

r 2 

rs 

function of r, s, and ct, 
7r 2 Jr 2 + (s 2 /cr 2 )] 
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k thermal conductivity of fluid 
H fluid viscosity 

p fluid density 

a channel aspect ratio 

Subscripts: 

B broad wall 

b bulk or mixed-mean value 
d fully developed 

e entrance value 


Q insulated wall or zero wall tempera- 
ture, internal heat generation 

q wall heat flux, no internal generation 

s sidewall 

T uniform wall temperature, no internal 
heat generation 

w value at wall 

Superscript: 

* entrance region 
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APPENDIX B 


TEMPERATURE DISTRIBUTION FOR INTERNAL HEAT GENERATION 
WITH WALL HEAT TRANSFER 

The steady- state temperature distribution for slug flow in a rectangular channel with 
internal heat generation and wall heat transfer is the solution to equation (2) subject to 
the boundary conditions expressed by equations (3). 

The problem has been examined from the standpoint of superposition of basic solu- 
tions in an earlier paper (ref. 28). This approach is a standard in heat- transfer analysis 
and need not be repeated here. The basic solutions may be found in reference 28; the re- 
sults have been used in the present report. 

This same problem can be treated by using the fact that the steady- state slug- flow 

convection problem is the same as the problem of unsteady- state heat conduction in a 

solid with the same cross-sectional shape and the same boundary conditions. Therefore, 

the results presented in reference 28, for the situation where a nongenerating liquid 

metal entering at a temperature t g flows through a channel with heat transfer q^ and 

q at the walls, might also have been derived from appropriate superposition of the solu- 
s 

tion of the problem of one- dimensional heat flow in a slab with uniform heat flux at the 
surface (ref. 29, p. 113). (This fact was called to the author’s attention by Mr. Aaron 
J. Friedland in a written discussion of reference 28 that was submitted to the 1966 Heat 
Transfer and Fluid Mechanics Institute. ) 

The temperature t_ (ref. 28) may be expressed as the sum t„ + t„ „ + t„ „, where 
tq g corresponds to the case of no internal heat generation, zero inlet temperature, 
heat flux q^ on the broad walls, and insulated sidewalls. The temperature t g cor- 
responds to the case of no internal heat generation, zero inlet temperature, heat flux q g 
on the sidewalls, and insulated broad walls. The problem is then reduced to the case of 
one -dimensional heat flow in a slab with uniform heat flux at the surface, and the general 
solution is obtained by adding the results for the two directions. The details are omitted, 
since the results present nothing new. 

Therefore, the solutions previously obtained in reference 28 for steady- state slug- 
flow convection could be found from consideration of the conduction of heat in the unsteady 
state. The advantage of the former approach to the problem, however-, lies in the fact 
that it is possible, in principle, to apply it to more difficult cases (e. g. , more realistic 
turbulent velocity profile; inclusion of turbulent diffusion as an additional mechanism for 
heat transfer) which are not easily dealt with by analogy to the unsteady- state heat- 
conduction problem. 
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APPENDIX C 


TEMPERATURE DISTRIBUTION FOR INTERNAL HEAT GENERATION 
WITH UNIFORM, EQUAL WALL TEMPERATURES 


The temperature t may be expressed as a sum tQ + t T , or in dimensionless form, 


t(x, y, z) - t, 


w 


*Q , *t ~ 


W 


*e ” *w 


't a - t t - t 
e w e w 


= T Q + T 1 


(Cl) 


The governing equations and boundary conditions for t q and Trp become 

OT Q ■ a \ , , Qa 2 A 

9x1 9y' 2 Sz' 2 ^e " V 


o 

11 

at y* = 0 and y T = 1 



t q = ° 

at z’ = 0 and z' = a 

► 

(C3) 

o 

11 

1 ? 

at x' = 0 

y 



3T t 

0X T 

3 2 T t 3 2 T t 

: H 

_ ,2 _ ,2 

9y' 3z' 


(C4) 

o 

ii 

H 

Eh 

at y T = 0 and y* = 1 



T t = 0 

at z’ = 0 and z f = a 

► 

(C5) 

T t = 1 

>1 

o 

II 

K 




The problems involving the temperatures Tq and Trj, have a definite physical 
meaning. The former is the case of a heat-generating liquid metal at t = 0 that enters 
a channel whose wall temperatures are also t = 0. The latter is the case of a nongener- 
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ating liquid metal at t that enters a channel whose wall temperatures are t w< The two 
problems are treated separately, and the results are combined to yield information for 
the general situation. 


Internal Heat Generation with Wall Temperatures t w = 0 


In the fully developed region the fluid flows isothermally (3Tq dx' = 0) and there- 
fore from equation (C2) T Q,d satisfies 


^Q.d , a %,d _ Qa 2 A 


3y 


,2 


3z* 


(t - t ) 
' e w' 


(C6) 


with the boundary conditions 


T Q,d=° 

at y' = 0 and y' = 1 

(C7) 

T Q,d=° 

at z’ = 0 and z' = a 

(C8) 


The boundary value problem posed by equations (C6) to (C8) is analogous to the prob- 
lem of torsion in rectangular bars (ref. 30). The analogy has been used also to treat the 
problem of laminar flow in a rectangular conduit (ref. 31). Thus, from the results given 
in reference 30, the temperature distribution in a heat- generating slug flow in a rectan- 
gular channel with wall temperatures t = 0 is given by 


T Q,dfr’’ z ’) 


oo 



n=l, 3, . . . 


(C9) 


The solution as given by equation (C9) has the advantage of avoiding the use of a double 
Fourier series for the temperature distribution. In principle, an extension of the tech- 
nique employed in reference 30 to obtain a single series solution of equation (C6) for solv- 
ing the analogous equation for torsion in rectangular beams makes possible a single 
Fourier series solution for the entrance region temperature Tq(x’, y', z’). The analytical 
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difficulties associated with satisfying the thermal entrance condition, however, lead to 
the mathematically intractable problem now presented. 

To determine the solution in the thermal entrance region, an entrance temperature 
TQ(x',y f ,z') is introduced such that 


T Q' T Q,d +T Q 


(CIO) 


After equation (CIO) is introduced into equation (C2) and it is taken into account that T Q,d 
satisfies equation (C6), the function Tq must satisfy the differential equation 


8T q - fzk , 3 ' t q 

9x1 3y’ 2 3z’ 2 


(Cll) 


with the boundary conditions 


1q-» 


t q=° 


at y’ = 0 and y’ = 1 


at z' = 0 and z’ = o 

j 


(C 12) 


At the channel entrance (x' = 0), the condition is 


T q (0, y*, z’) = 0 = T q> d (y', z’) + Tq( 0, y',z’) 


or, by rearranging, 


Tq( 0, y f , z') = 


Qa 
4a^ k 


3 t 0 
7T e 


w 



n 


cosh — 
1 


r-i 

V 2 


cosh — 
2 ct 


sm 


n7rz* 


n=l, 3, . . 


(Cl 3) 


Equation (C13) suggests taking T^(x*, y T ,z') in the form of a series: 
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(Cl 4) 



n=l, 3, . . . 


_j3^x’ 

A sin w e n 
n _ n 


where Aj, Ag, . . . are constant coefficients and Wp Wg, . . . are functions of y* 
only. Substituting equation (C14) into equation (Cll) and taking into account the boundary 
conditions (C12) yield the general expression for the entrance temperature: 


n „2 go 

Qa \ \ o 

\ » -px' 

Ti(x',y',z')= - — \ A_ sin nn-y’ sin e n (C15) 

Q *e-‘w L, 

n=l, 3, . . . 

2 2 2 2 

where = n 7 t (1 + 1/cr ). The coefficients A n have yet to be determined. Applying 
the boundary condition at the channel entrance as given by equation (C13) results in 



n=l, 3, . . . n=l,3,... 


(C16) 

The solution of equation (C16) for the constant coefficients A n does not appear possible, 
and it is concluded that the formulation of the temperature distribution Tq in terms of 
a single Fourier series does not appear to offer any significant advantages. The temper- 
ature distribution, on the other hand, is susceptible to a straightforward analysis in 
terms of a double Fourier series. Therefore, the problem is now treated from the stand- 
point of double Fourier series. The solution of equation (C6) can be expressed as a 
double sine series: 


t q>’’ z, > = 


Qa 


K(t e 




oO oo 

zz 

m=l n=l 


z * 

C mn sin(m7ry’)sin nn — 

a 


(C17) 
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The individual terms and the entire series vanish on the boundaries. Hence, the boundary 
conditions (eqs. (C7) and (C8)) are satisfied. 

The next task is to determine the Fourier expansion coefficients C mn so that equa- 
tion (C17) will satisfy the governing equation (C6). Substituting equation (C17) into equa- 
tion (C6), multiplying by a typical product of sines, and finally integrating over the cross 
section of the channel give 


C 


mn 


rnnjr 


16 


2p2 


mn 


(C18) 


2 2T 2 2 21 

where r mn = n m + (n /cr ) and where m and n are odd integers. If either m or 
n or both are even numbers, C mn = 0. 

The entrance temperature Tq must now be determined, and it must satisfy equa- 
tions (Cll) and (C12). The boundary condition corresponding to the requirement that the 
fluid enter the channel with a temperature value t = 0 is now 


2 oo oo 

Tq( 0, y', z') = -T q d (y*, z») = — Qa V V C mn sin(m ff y0sinfn7r^ (C19) 

’ K ' t e ‘ V L t-n 1 \ °t 

m=l n=l 

Equation (Cll) and the associated boundary conditions can be solved by the 
separation- of- variables technique. Let 

t q = <Pi( x ')^2(y')<P 3 (z T ) (C20) 

where cp^x'), <P 2 (y')> and <p 3 (z*) are functions of x% y T , and z ' , respectively. Substi- 
tution of equation (C20) into the differential equation (Cll) and consideration of the bound- 
ary conditions (eq. (C12)) give the solution for Tq in the form of a double sine series: 


°© °° ^ 

t « = Z Z Drs sin(r ' y,)sin f' 7) e ’ rrsX ' (c2i> 

r=l s=l 

where r^ g = 7r 2 |r 2 + (s 2 /cr 2 )j . The coefficients D rg remain to be determined. Applying 
the boundary condition at the channel entrance as given by equation (C19) yields 
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m=l n=l 


mn 


sin(m3ry')sin 



— \ ' \ ' D „ sin(rjry')sin/s7r 

- 1 ) Z, rs 
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r=l s=l 


V 


3 


(C22) 


The equality of the two doubly infinite series in equation (C22) is interpreted to be 

equivalent to the results r = m, s = n, and D = -C . Hence, 

’ rs m ’ 



Qa 


«(t c 


OO QQ 

bEE 

m=l, 3 n=l, 3 


C mn sin(m7ry')sin|raT — ]e 
' or 


-r 2 x» 
mn 


(C23) 


Combining the results of the preceding sections in accordance with equation (CIO) 
gives an expression for the temperature distribution that applies in both the entrance and 
fully developed regions: 


T Q (x',y%z») 



16 (l 

m - 2r L ' 



sin(m7ry’)sin n?r — 
cr 


(C24) 


The local bulk temperature 


L Q,b 


of the liquid metal is given by 



(C25) 


Introducing equation (C24) into equation (C25), integrating and then rearranging yield 
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oo oo 



m=l, 3, . . . n=l, 3, . . . 


(C26) 


Wall Temperatures t w Without Internal Heat Generation 

The dimensionless temperature T T (x',y',z') is the solution to equations (C4) and 
(C5). (The problem is treated in ref. 12 from the standpoint of a Fourier analysis.) The 
temperature field is described by the double infinite Fourier series (ref. 12) 

o° °o 2 

Z - 16 sin(m77y')sin (m — ) e mn (C27) 

mmr^ ' 

m=l n=l 

where m = 1, 3, 5, . . . and n = 1, 3, 5, . . . . The local bulk temperature T T b 
may be found as follows: 

^ = T t , = 1 f° f X T t dy’ dz’ (C28) 

t e -t w T ’ b a JO Jo T 

Substituting the temperature distribution from equation (C27) into equation (28) and inte- 
grating result in 



OO OO 



m=l n=l 

Combining the solutions for t q and Trp yields the general result 


(C29) 
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